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摘要 :干旱 指数 的 区 域 适用 性 是 准确 表征 区 域 干旱 的 重要 前 提 , 本 文 以 中 国 东北 地 区 为 典型 研究 区 ,探讨 标准 化 降水 蒸 散 指数 
(SPEI) 在 该 地 区 应 用 的 有 效 性 。 基 于 研究 区 90 个 气象 台 站 的 逐日 气象 资料 ,计算 1961 一 2014 年 多 时 间 司 广 的 SPEI 指数 。 从 
Kolmogorov-Smirnov( KK-S) 拟 合 优 度 检验 SPEI 与 典型 干旱 事件 核准 、SPEI 与 农作物 受 旱灾 面积 及 与 土壤 湿度 相关 性 分 析 等 方 
面 ,验证 SPEI 指数 在 东北 地 区 的 适用 性 。 分 析 结 果 表明 :1) 东北 地 区 多 时 间 尺 度 的 累积 玲 分 方 缺 量 符 合 Log-logistic 分 布 ， 
SPEI 指数 在 东北 地 区 的 应 用 具备 数学 统计 理论 基础 ;2) 生长 季 平 均 SPEI 值 与 黑龙 江 省 .吉林 省 和 辽 字 省 农作物 受 旱 灾 面 积 比 
例 均 呈 极 显 著 负 相关 (P<0.01) ;3) 在 1.3.6 和 12 个 月 尺度 下 ,SPEI 与 土壤 湿 语 呈 显 著 正 相关 (<0.05) 的 站 点 比例 分 别 为 
90.2% .92.16% .90.2% 和 88.24%。 综 上 所 述 ,SPEI 指数 不 仅 满 足 数 学 理论 统计 的 要 下 而 且 与 干旱 灾情 数据 和 土壤 水 分 监测 值 
均 具有 极度 的 关联 性 ,说 明 其 在 东北 地 区 干旱 预测 和 定量 化 研究 中 具有 较 好 的 适用 性 。 
关键 词 :标准 化 降水 蒸 散 指数 ;K-S 检验 ;土壤 湿度 ;干旱 受灾 面积 ;干旱 事件 


= 


全 


Applicability analysis of SPEI for drought research in Northeast China 
SHEN Guoqiang" “，ZHENG Haifeng * , LEI Zhefifeng 


1 Northeast Institute of Geography and Agroecology, Chinese Acdderiy of WSciences, Changchun 130102, China 
2 University of Chinese Academy of Sciences, Beijing 100049, China 


Abstract: Determining the applicability of a drought index to a specific area is a prerequisite for drought quantification and 
mapping on a regional scale. In this Wpaper, we explored the effectiveness of the Standardized Precipitation 
Evapotranspiration Index ( SPEI)Sin northeast China. SPEI values at a time scale of 1, 3, 6, and 12 months were calculated 
using a daily climate dataset thatlhad been collected by 90 meteorological stations in northeast China from 1961 to 2014. The 
climate variables analyzed included daily temperature, precipitation, atmospheric pressure, wind speed, relative humidity, 
and sunshine duration. The applicability of the SPEI for drought quantification in northeast China was investigated using 
Kolmogorov-Smirnov (K-S)》 test, validation analysis from specific drought events and soil moisture, and correlation analyses 
between, SPED and the statistics of the area suffering from drought.The results showed that our null hypothesis was rejected 
by the 下 -SiTest atfa 0.05 significance level for each station, and that the Log-logistic distribution matched very well to the 
water stirplus or deficit series for all four time scales at each station. The drought-affected areas showed a significant negative 
relationship to the SPEI during the growing season ( April to September) at these four time scales during 1971 一 2013. The 
index of SPEIO1 was found to be more suitable for predicting the time, location, and intensity of drought events. There was 
a strong positive correlation between soil moisture at a depth of 20 cm and the SPEI for most of the climate stations ( 90.2% 


of stations for SPEIO1, 92.16% of stations for SPEI03, 90.2% of stations for SPEI06 and 88.24% of stations for SPEI12 ) . 
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In addition, the number of climate stations decreased as the time scale increased. Our results indicate that the SPEI is 
suitable for quantifying droughts, and therefore, should be widely used in scientific drought prediction research in Northeast 


China. 
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干旱 指数 是 定量 表征 干旱 事件 的 重要 指标 ,常用 的 干旱 指数 有 降水 距 平 百分率 、 土 壤 湿 度 、 燕 发 量 z 降 水 
量 、Palmer 干旱 指数 (Palmer drought severity index, PDSI) 和 标准 化 降水 指数 (Standardized precipitationindex， 
SPD) 等 … ,其 中 PDSI 和 SPI 应 用 较为 广泛 。PDSI 指数 的 优点 在 于 考虑 了 温度 和 前 期 天 气 条件 对 干旱 的 影 
响 ,非常 适用 于 干旱 对 全 球 变 暖 响应 的 研究 ,但 其 参数 获取 困难 且 计 算 复杂 ,加 之 时 间 尺 谋 固 害 @9 一 12 个 
月 之 间 ) ,因此 无 法 有 效 地 应 用 于 干旱 的 多 时 间 尺 度 研 究 中 中 。SPI 指数 虽然 适用 年 多 时 间 尺 度 研 究 ,但 其 仅 
以 降水 量 作为 干旱 定量 的 唯一 参数 ,忽略 了 温度 对 干旱 形成 的 重要 作用 “| /标准 站 降 水 蒸 散 指数 
( Standardized precipitation evapotranspiration index ,SPEI) 不 仅 充 分 考虑 了 气温 对 干旱 的 影响 ,而 且 综 合 考 
虚 了 干旱 的 多 时 间 尺 度 。SPEI 自 2010 年 被 提出 以 来 ,被 广泛 应 用 于 定量 琅 量 的 研究 中 ,例如 SPEI 与 气候 因 
子 的 相互 关系 研究 !'9 不同 干旱 指数 的 对 比分 析 ' ”区域 干 时 时空 分 布 特征 分 析 '" 六 散 量 算法 差异 对 SPEI 
的 影响 中 以 及 SPEI 的 多 时 间 尺 度 研究 '] 。 鉴 于 干旱 的 复杂 性 广泛 人 性 和 时 空 差 异性 ,干旱 指数 的 区 域 适 用 
性 成 为 干旱 研究 的 前 提 , 如果 不 根 据 区 域 的 实际 情况 进行 计算 ,而 将 SREI 应 用 于 定量 干旱 事件 ,容易 导致 结 
论 偏差 "”。 然 而 ,在 相关 SPEI 的 研究 中 , 仅 有 少数 学 者 关注 SPET 指 数 的 前 提 和 适用 性 。 这 些 研 究 大 多 以 年 
为 时 间 尺 度 或 以 单个 月 份 为 研究 尺度 (例如 7 月份 和 12 月 份 ) 探 讨 SPEI 在 全 球 或 全 国 尺 度 的 适用 性 [5 ， 
并 没有 考虑 SPEI 指数 的 假设 和 前 提 。 由 于 区 域 干旱 精确 定量 化 的 需求 在 于 一 干旱 指数 要 准确 反映 任 一 月 
份 、 任 一 站 点 的 水 分 收文 状况 ,因而 ,人 迫切 需要 以 月 为 时 间 斥 度 ,验证 该 指数 在 区 域 尺 度 上 所 有 站 点 的 适用 性 。 

自 20 址 纪 90 年 代 中 期 以 来 ,东北 地 区 干旱 化 趋势 不 断 增强 "9 ,生态 系统 的 干旱 脆弱 性 随 之 加 大 ,农业 
发 展 和 生态 系统 平衡 面临 巨大 威胁 。 精确 定量 东北 地 区 的 和 干旱 ,对 该 区 的 干旱 适应 性 研究 与 管理 具有 
重要 意义 。 到 目前 为 止 ,SPEI 指数 在 定量 东北 地 区 干旱 方面 的 有 效 性 和 适用 性 尚 不 明确 。 本 文 将 从 SPEI 的 
理论 假设 和 实际 定量 验证 入 手 ,应 用 正 -$ 检验 和 相关 性 分 析 等 方法 ,探讨 (1) SPEI 指数 在 东北 地 区 的 数学 统 
计 理 论 基 础 ;(2) 干 旱 事件 .旱灾 面积 和 土壤 湿度 对 SPEI 指数 的 验证 效果 。 


1 研究 区 与 数据 


1.1 研究 区 概况 

研究 区 为 中 国 东 北上 籽 区 (38"43'N 一 53"33'N，118"50' 上 一 135"2E) (图 1) ,行政 上 包括 黑龙 江 省 .吉林 省 
和 辽宁 省 ,总 面积 约 84.53 万 km 。 该 区 面积 广阔 ,地 形 以 山地 和 平原 为 主 , 其 中 ,东北 平原 黑土 资源 丰富 , 土 
壤 肥 沃 ,是 我 国 莉 要 的 粮食 基地 。 该 区 主要 气候 类 型 为 温带 季风 气候 ,年 平均 气温 为 5.4%C ,由 南 向 北 依次 跨 
越 暧 温带 、 申 温带 和 寒 温带 ;年 平均 降水 量 约 为 600mm ,由 东 向 西 逐 渐 减 少 ,东部 为 湿润 区 ,西部 为 半 湿 润 区 。 
研究 区 干旱 频 发 ,尤其 是 进入 21 世纪 之 后 ,干旱 事件 发 生 更 加 频繁 .持续 时 间 更 长 '”" ,这 严重 制约 东北 地 区 
农业 生产 的 发 展 。 据 农业 部 种 植 业 管理 司 统计 ,1971 一 2013 年 ,东北 三 省 平均 每 年 受 旱 灾 面 积 达 4x104kms 
(上 由 地 区 播种 总 面积 的 22.4% ) ,其 中 黑龙 江 省 .吉林 省 和 辽宁 省 受灾 面积 分 别 为 1.87x10' km” 1.14x10 km 和 
9.9x10 km ,分别 占 各 自 农 作物 播种 总 面积 的 20.1% .24.9% 和 25% 。 
1.2 ”数据 来 源 

本 文 所 用 气象 数据 来 源 于 中 国 气象 局 发 布 的 中 国 地 面 气 候 资 料 日 值 数据 集 (V3.0) ,包括 研究 区 90 个 气 
象 站 点 (黑龙 江 .吉林 和 辽宁 分 别 占 31.27 和 32 个 )1961 年 1 月 1 日 到 2014 年 4 月 30 日 的 逐日 气象 要 素 资 
料 ( 气 温 、 降 水 量 、 气 压 、 风 速 .相对 湿度 和 日 照 时 数 )。 这 些 气 象 数 据 经 过 中 国 气 象 局 严格 地 质量 控制 ,数据 
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的 实 有 率 和 正确 率 在 99% 以 上 。 

典型 干旱 事件 资料 来 源 于 中 国 气 象 灾 害 年 鉴 。 本 
文 整理 了 2004 一 2010 年 东北 地 区 影响 较 大 的 8 次 干旱 。 
事件 的 时 间 .范围 和 强度 。 

农作物 受 旱 灾 面 积 资料 来 自 农 业 部 种 植 业 管理 司 
灾情 数据 库 (202.127.42.157/moazzys/zaiqing.aspx) ,本 
文 分 析 了 1971 一 2013 年 东北 各 省 的 农作物 播种 面积 以 
及 受 旱 灾 面 积 。 

土壤 湿度 数据 来 源 于 中 国 气象 数据 网 ( data. cma. 
cn) 发 布 的 中 国 农作物 生长 发 育 和 农田 土壤 湿度 旬 值 
数据 集 , 数 据 的 时 间 范 围 为 1991 一 2013 年 。 为 了 保证 
数据 的 完整 性 并 减少 灌溉 措施 对 土壤 湿度 的 影响 ,文中 。 
保留 了 51 个 非 水 田 站 点 的 土壤 湿度 (10 一 20cm 深 处 ) 
数据 ,用 算数 平均 法 将 旬 值 转 为 月 值 。 
2 研究 方法 
2.1 SPEI 指数 计算 

SPEI 以 月 水 分 亏 缺 量 ( 即 降水 量 减 去 蒸 散 量 ) 来 表 0 ，150 300km 
示 水 分 的 僵 余 或 赤字 。 基 于 “历史 同月 的 累积 水 分 亏 
缺 量 服从 Log-logistic 分 布 ” 这 一 假设 ,计算 分 布 概率 密 


120° 123° 130° 133° 140°E 


度 函 数 和 累计 概率 ,然后 转化 成 标准 正 态 分 布 进而 求 a 
得 ! o 具体 计算 步骤 如 下 : Fig.l Study area and Meteorological stations 


第 一 步 计算 潜在 蒸 散 量 。 本 文采 用 FAO Penman- 
Monteith 法 ' 站 计算 潜在 蒸 散 量 ,公式 如 下 } 


0.408A(R ~ G) +y (e, 一 eu) 


900 
T4273 
A+y(l +0.34n,) 
式 中 ,PET 为 潜在 蒸 散 量 (Cmigid) , RR, 为 地 表 净 辐射 (MJ md "1!) ,6 为 土壤 热 通 量 (MJ md!') ,了 为 日 平均 
气温 (% ) , A 为 饱和 水 洽 压 曲线 斜率 ( kPa/%C ) , y 为 干 湿 表 常数 (kPa/%C ) ,为 2m 高 处 的 风速 (m/s)，, e， 
为 饱和 水 汽 压 (kBa) ,e ,为 实际 水 汽 压 (kPa) 。 
第 二 步 计算 逐 月 降水 与 蒸 散 的 差 值 , 即 水 分 亏 缺 量 : 
D,=P, - PET 
式 趾 ,D, 为 水 分 纪 缺 量 ,P, 为 月 降水 量 ,PET, 为 月 潜在 燕 散 量 ; 并 根据 线性 递减 权重 ' "方案 建立 不 同时 间 尺 度 
的 办 积 水 分 亏 缺 量 序列 . 


PET = 


:s_ [2(i+1) 
DE |p, PET,) | n>1 
式 中 ,为 时 间 尺 度 ( 月 ) ,n 为 计算 次 数 ; 
第 三 步 采用 Log-logistic 概率 密度 函数 拟 合 所 建立 的 水 分 亏 缺 量 序列 : 


A -人 a SDT 


式 中 ,a 为 尺度 参数 ,8 为 形状 参数 ,y 为 位 置 参数 ,可 通过 线性 矩 的 方法 拟 合 获得 。 累 积 概率 可 通过 分 布 函数 
计算 : 
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第 四 步 对 拟 合 的 水 分 亏 缺 量 序列 进行 标准 正 态 分 布 转换 ,获得 对 应 的 SPEI: 
spEI 2W- I ,WW= /2In (Pp) 
1l+dW+d, W +d,W 

当 P<0.5 时 ,P=1 -F(x); 当 P>0.5 时 ,P=1 -PP, 同时 SPEI 的 符号 逆转 , 式 中 其 他 常数 项 分 别 为 Gy= 
2.515517, C, =0.802853, C, = 0.010328 ，d| = 1.432788 ，d, = 0.189269，d; = 0.001308 。 蒸 散 量 和 SPEI 值 均 R 语 
言 SPEI package 中 进行 ,本 文 计 算 了 1.3.6 和 12 个 月 尺度 的 SPEI ,分 别 记 为 SPEIO1 .SPEI03 SPEI06 和 SPEI12。 
2.2 K-S 拟 合 优 度 检验 

K-S 检验 '” 是 一 种 拟 合 优 度 检验 方法 ,利用 样本 数据 推断 样本 来 自 的 总 体 是 否 与 菜 一 理论 分 布 有 显著 
差异 。 其 基本 思路 是 :用 (x) 表示 各 样本 观测 值 在 理论 分 布 上 出 现 的 理论 累积 概率 值 , S(x) 表 示 各 样本 观 
测 值 的 实际 累计 概率 值 , 计 算 检验 统计 量 . 

D=max| S(x) —- F(x) | 

当 零 假设 成 立时 ,D 统计 量 服从 Kolmogorov 分 布 ,计算 对 应 的 概率 p 值 s 若 了 加 于 显著 性 水 平 a, 则 应 拒 
绝 零 假 设 , 认 为 样本 来 自 的 总 体 与 指定 的 分 布 有 显著 差异 ; 否则 ,不 能 拒绝 零 假 设 , 记 为 二 者 无 显著 差异 。 为 
了 验证 东北 地 区 累积 水 分 亏 缺 量 是 否 符合 Log-logistic 分 布 , 从 而 判断 SPEI 在 东北 地 区 的 应 用 是 否 具 备 数理 
统计 基础 ,本 文 对 不 同时 间 尺 度 下 的 累积 水 分 亏 缺 量 序列 与 Log-logistia 分 布 进 行 了 K-S 检验 。 
2.3 ”相关 性 分 析 

为 了 验证 SPEI 在 表征 农作物 受 旱灾 方面 的 适用 性 ,本 闫 分 别 分 析 了 1971 一 2013 年 黑龙 江 省 .吉林 省 和 
辽宁 省 生长 季 (4 一 9 月 ) 平 均 SPEI 值 与 当年 作物 受 旱灾 面积 比例 的 相关 性 (各 省 的 样本 量 均 为 43a, 当 n=43 
时 ,rs-oo 的 临界 值 为 0.389) 。 

同时 ,为 了 验证 SPEI 在 反映 土壤 干 湿 状 狐 方 面 的 适用 性 ,本 文 以 站 点 为 检验 单元 ,分 析 了 土壤 湿度 与 
SPEI 指数 的 相关 性 (各 站 点 的 样本 量 均 为 134 个 月 , 当 n=134 时 ,ro 的 临界 值 为 0.17,r -ou 的 临界 值 为 
0.222)。 文 中 相关 系数 7 均 指 Pearson4 简 单 相关 系数 。 


3 结果 分 析 


3.1 基于 K-S 检验 的 SPIE 指数 理论 验证 

累积 水 分 亏 缺 量 序列 与 Logslogistic 分 布 的 K-S 检验 结果 (图 2) 表 明 , 在 不 同时 间 尺 度 下 ,K-S 检验 的 概 
率 p 值 均 大 于 显著 性 水 平 a=0.05 ,没有 足够 理由 拒绝 零 假 设 ( 零 假 设 为 累积 水 分 气 缺 量 样本 来 自 Log-logistic 
分 布 总 体 )。 因 此 ,可 以 判定 东北 地 区 累积 水 分 亏 缺 量 序列 与 Log-logistic 分 布 无 显著 差异 , 即 该 序列 来 自 
Log-logistic 分 布 总 体 。 在 不 同 尺 度 不同 站 点 和 不 同月 份 条 件 下 ,Log-logistic 分 布 函数 能 均 较 好 地 拟 合 累 积 
水 分 气 缺 量 序列 * 由 此 表明 ,基于 该 分 布 的 SPEI 指数 在 东北 地 区 干旱 表征 上 具备 数学 统计 理论 基础 。 
3.2“ 典 型 干旱 事件 对 SPEI 指数 的 核准 分 析 

典型 后 旱 事件 与 同期 SPEI 核准 的 结果 ( 表 1) 表 明 ,SPEI 与 干旱 事件 的 时 间 、 地 点 和 强度 吻合 度 较 高 。 
例如 ,在 气象 灾害 年 鉴 中 的 记录 为 “2006 年 4 一 5 月 黑龙 江 齐 齐 哈 尔 ` 大庆 和 绥化 等 地 旱情 严重 (中 到 重 
度 ) ”, 经 查 ,齐齐哈尔 .大庆 绥化 和 海伦 2006 年 5 月 的 SPEI 值 分 别 为 -1.3,-1,-1.2 和 -1.7, 达 到 中 度 甚至 
重度 干旱 ,与 年 鉴 记载 相 符 。 另 据 记 载 ,2009 年 6 月 21 日 至 11 月 上 名, 辽宁 和 吉林 出 现 了 大 范围 的 中 到 重度 
气象 干旱 ,吉林 农作物 受 旱 面 积 2.85x10‘km? , 重 旱 面 积 1.32x10*km? ,8.2 万 农村 人 口 .6.4 万 头 大 牲畜 因 早 发 
生 饮 水 困难 ;其 中 ,长 岭 县 普遍 遭受 伏旱 ,绝收 面积 达到 160 km? ,经 查 ,长 岭 站 2009 年 7 一 9 月 的 SPEI 分 别 为 
-0.8 -1.8 和 -1.6, 其 中 8.9 月 达到 重度 干旱 。 由 此 可 见 ，SPEI 能 较 好 地 反映 东北 地 区 干旱 事件 的 时 间 、 地 
点 及 强度 。 
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图 2 K-S 检验 概率 p 值 统计 图 ( 虚线 为 显著 性 水 于 ga=0,05) 
Fig.2 Thep values of K-S Test for water surplus or deficit ( dashedjline was the significance level: «=0.05) 


SPEI: 标准 化 降水 蒸 散 指数 Standardized precipitatiom evapotranspiration index 


表 1 典型 干旱 事件 与 同期 (SPETI 核准 (2004 一 2010 年 ) 
Table 1 Verification of the DroughtlEvents and SPEI in 2004 一 2010 


持续 时 间 范围 早 灾 程度 SPEI 
Duration Site Hazard intensity 
i 黑龙 江西 部 昌 图 5 月 : -1.3; 
AP 2004 年 3 月 一 6 月 中 旬 吉林 西部 重 旱 ,部 分 特 旱 富裕 6 月 : -1.6; 
辽宁 北部 白城 3 一 6 月: -1.4, -1.1,-0.7, -1.4; 
此 A 和 齐齐哈尔 6 一 7 月 ; -2, -1.1; 
六 年 5 月 黑龙 浅 西 部 重 早 SS a 。 
一 2004 年 5 月 下 旬 一 7 月 上 旬 黑龙 汪 西 音 重 早 龙 江 6_7 月 : -19 -18， 
齐齐哈尔 5 月 : -1.3; 
ee 大 庆 5 月: -1; 
年 4 月 一 丑态 T 中 两 癌 中 到 重度 于 里 y 
2006 年 4 月 一 5 月 黑 郑 江 中 西部 到 重度 干旱 绥化 5 月 ; _1.2， 
海伦 5 月; -1.7; 


宝 清 7 月 : -1.2; 


2007 年 6 月 11 日 -8 月 10 日 三 江平 原 旱情 严重 富 锦 6 一 7 月 : -1.1, -1.6; 
佳木斯 6 一 7 月 : -1.2, -2.2; 
笠 玛 1 一 2 月 : -1.5，-1.33; 
2008 年 1 用 3 月 中 旬 东北 全 境 1951 年 以 来 降水 最 少 1 月 东北 平均 : -0.62 
2 月 东北 平均 : -0.75 
上 寺 东 9 一 10 月 : -1.6,-1.3; 
黑龙 江 人 
2008 年 9' 月 下 旬 一 10 月 下 旬 洒 守 降水 偏 少 30% 一 90% 桓 仁 9 一 10 月 : -2,-1; 
齐齐哈尔 9 月 ; -1 2 


四 平 7 一 9 月: -1.3,-1.2,-1.5; 
农安 7 一 9 月 : -0.8,-1.5,-1; 


四 辽宁 1952 年 以 来 最 严重 长 岭 7 一 9 月 : -0.8,-1.8,-1.6; 
ee 夏秋 连 早 开 原 7_9 月 ; -1,-1.6,-18; 


朝阳 7 一 9 月 : -1.6,-1.7,-1.3; 
义 县 7 一 9 月 :-1.4,-1.4, -1.6; 
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持续 时 间 范围 旱灾 程度 SPEI 
Duration Site Hazard intensity 
大 庆 9 月 H =133 
肇 州 9 月 : —1.4; 
秘 源 9 月 : -1.8; 
六 : ; 
2010 年 9 月 ne 中 到 重度 干旱 哈尔滨 9 月 : -1.6; 
可 六 
黑 乾 安 9 月， -1.4; 
通 榆 9 月 过 一 1.4; 
大 安 9 月 : -1.8 


SPFI: 标准 化 降水 蒸 散 指数 Standardized precipitation evapotranspiration index 


3.3 干旱 受灾 面积 对 SPEI 指数 的 验证 分 析 

农作物 受 旱灾 面积 比例 与 生长 季 平 均 SPEI 值 的 相关 性 分 析 结 果 显 示 , 各 省 农作物 受 旱灾 面积 比例 与 
SPEIO1 SPEI03 .SPEI06 和 SPEI12 的 相关 系数 在 -0.52 与 -0.76 之 间 ,SPEI 指数 与 受 早 帘 面 积 呈 极 显 著 负 相 
关 关 系 (P<0.01 ) ,表明 干旱 越 严 重 ,SPEI 值 越 小 , 受 旱灾 面积 越 大 。 同 时 ,研究 发 现 , 在 相同 的 时 间 尺 度 下 ， 
受 早 灾 面 积 比例 与 SPEI 的 相关 性 在 各 个 省 份 之 间 也 存在 差异 ,相关 性 最 强 的 是 辽宁 省 $ 其 次 是 吉林 省 ,黑龙 
江 省 相关 性 最 弱 。 


ee 生长 季 平 均 SPEI 一 一 一 一 一 一 受灾 面积 
黑龙 江 省 吉林 省 辽宁 省 


r=-0.76 (P<0.01) 


SPEIO1 


r= -0.67 (P<001) 


’ 
由 


SPEI03 


标准 化 值 Standardized value 
1 


SPEI06 


SPEI12 


1970 1980 1990 2000 2010 1970 1980 1990 2000 2010 1970 1980 1990 2000 2010 
年 份 Year 


图 3 东北 三 省 农作物 受灾 面积 比例 与 生长 季 平均 SPEI 对 比 


Fig.3 The comparison between drought crop area ratio and the mean SPEI 


3.4 土壤 湿度 对 SPEI 指数 的 验证 分 析 

SPEI 指数 与 土壤 湿度 的 相关 性 分 析 结 果 ( 图 4) 表 明 , 在 1.3.6 和 12 个 月 尺度 下 ,51 个 研究 站 点 中 分 别 
有 46 .47.46 和 45 个 站 点 通过 了 相关 系数 的 显著 性 检验 (P<0.05) ,通过 率 分 别 达 90.2% .92.16% .90.2% 和 
88.24% ;其 中 ,分 别 有 46 .45 .44 和 40 个 站 点 呈现 极 显 著 相 关 (P<0.01) ,通过 率 分 别 为 90.2% .88.24% 、86.27 
和 78.43%。 同 时 ,研究 发 现 站 点 显著 性 检验 的 通过 率 随 着 SPEI 时 间 尺 度 的 增 大 而 降低 。 通 过 了 显著 性 检验 
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(P<0.05) 的 所 有 研究 站 点 ,其 土壤 湿度 与 SPEI 的 相关 系数 取 值 范围 为 0.26 一 0.74 ,平均 值 为 0.47, 表 明 土 壤 
湿度 与 SPEI 呈 显 著 正 相关 ,SPEI 值 越 大 ,土壤 湿度 越 大 。 


201706.00460v1 
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4 SPEI 与 土壤 湿度 相关 分 析 空 间 分 布 


Fig.4 Spatial distribution of correlation analysis of SPEI and soil moisture 


人 ED) 根 据 SPEI 指数 的 定义 ,SPEI 反映 的 是 某 月 累积 水 分 气 缺 量 在 历史 同期 累积 水 分 亏 缺 量 序列 中 的 分 
布 概率 ,所 以 应 当 用 Log-logistic 分 布 对 1 一 12 月 的 累积 水 分 气 缺 量 序 列 进行 逐 月 地 拟 合 ; 而 Vicente- 
Serranol3l 和 庄 少 伟 : 呈 在 验证 SPEI 的 适用 性 时 ,是 对 整个 研究 时 期 每 年 的 水 分 亏 缺 量 序列 进行 拟 合 检验 ,这 
种 以 年 为 时 间 尺 度 的 验证 方式 掩盖 了 不 同月 份 之 间 的 水 分 亏 缺 量 差 异 ,偏离 了 SPEI 的 前 提 假 设 , 据 此 得 出 的 
结论 有 待 商 检 。 王 林 !51 验 证 了 在 1 个 月 尺度 下 ,东北 地 区 7 月 份 和 12 月 份 的 水 分 亏 缺 量 符合 Log-logistic 分 
布 ,其 结果 与 本 文 相符 ,但 未 考虑 前 期 降雨 对 干旱 的 影响 及 多 时 间 尺 度 SPEI 计算 的 重要 性 。 为 了 确保 SPEI 
指数 能 精确 定量 地 区 干旱 ,本文 对 东北 地 区 各 个 站 点 的 水 分 亏 缺 量 序列 进行 逐 月 地 拟 合 ,并 逐一 进行 K-S 检 
验 , 证 明了 SPEI 指数 在 东北 地 区 具有 较 好 的 数学 统计 理论 基础 。 
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(2)SPEI 的 理论 适用 性 取决 于 所 选 的 理论 分 布 。Vicente-Serranol53l 提出 SPEI 指数 时 ,对 比 了 Pearson 
亚 .Lognormal 广义 极 值 分 布 ( generalized extreme value，GEV) 和 Log-logistic 四 种 理论 分 布 对 水 分 亏 缺 量 经 验 
分 布 的 拟 合 效果 ,最终 认 为 Log-logistic 更 符合 其 研究 站 点 的 水 分 收 支 实际 情况 。 王 澄海 汪 和 王 芝 兰 :入 对 中 
国 的 降水 量 序列 进行 了 GEYV 分 布 的 拟 合 ,并 构建 了 GEV 指数 ,虽然 忽视 了 东 散 对 干旱 的 影响 ,但 为 干旱 指数 
的 地 区 适用 性 研究 开拓 了 新 的 视野 , 即 应 当 根 据 研究 区 水 分 收 支 特点 建立 区 域 适 用 的 分 布 模型 ,并 据 此 构建 
干旱 指数 。 此 外 ,由 于 不 同 检验 方法 功效 不 同 ,所 以 还 应 选取 合适 的 方法 来 检验 水 分 亏 缺 量 理论 分 布 与 经 验 
分 布 的 拟 合 优 度 , 如 K-S,，Anderson-Darling( A-D) 和 Shapiro-Wilk(S-W) 等 滞 。 综 上 ,在 验证 干旱 指数 的 区 域 
适用 性 时 ,有 必要 对 多 种 理论 分 布 进行 严格 的 拟 合 优 度 检验 ,确保 研究 结论 在 统计 学 上 的 合理 性 。 本 研究 所 
使 用 的 Log-logistic 分 布 是 干旱 研究 中 较为 常用 的 分 布 函数 ,并 且 已 被 证 明 其 与 东北 地 区 水 分 亏 缺 量 的 拟 合 优 
度 较 好 。 然 而 ,是 否 有 拟 合 效果 更 优 的 理论 分 布 ? 这 也 是 干旱 精确 定量 化 研究 的 重要 课题 之 会 4 

(3) 在 同一 时 间 尺 度 下 ,东北 三 省 受 旱 灾 面 积 与 


25.0 


SPFI 指数 的 相关 性 从 强 到 弱 依次 为 :辽宁 > 吉林 > 黑龙 。 x 225| 

江 。 出 现 这 个 结果 的 主要 原因 可 能 是 受 农业 种 植 结构 村 这 175 | 

变化 的 影响 。 陈 莉 529 认 为 水 稻 栽培 多 以 水 利 灌溉 为 基 车 号 150| 

而 ,受降 水 影响 不 大 ,因此 在 分 析 多 年 农作物 干旱 面积 县 ioo 上 

时 ,应 排除 水 稻 面 积 增加 的 影响 ,只 考虑 旱 作 农作物 的 。 < so[ 

播种 面 从 。 本 文 认为 ,水 稳 的 生理 抗旱 性 虽然 低 于 其 他 《 生 Le 1980 1985 1950 1995 2000 2005 2010 
作物 ,但 由 于 灌溉 条 件 成 熟 ,水 稻 受 干旱 的 影响 比 其 他 年 份 Year 


作物 小 , 即 同 样 的 干旱 强度 下 ,其 他 作物 比 水 稻 更 易 受 0 
灾 。 也 就 是 说 ,SPEI 指数 对 玉米 ` 小 麦 等 农作物 受 早 情 Fig.5 The variation of area ratio of paddy in the Northeast China 
况 反 映 灵敏 ,而 对 水 稳 受 旱情 况 反映 不 灵敏 。 由 图 5 可 
见 ,1971 一 2013 年 黑龙 江 省 水 称 面 积 由 1.65xJozkm 扩大 到 3.18x10km?, 占 农作物 总 面积 比例 由 2% 上升 到 
26% ,绝对 面积 和 所 占 比 例 在 2008 年 之 后 鬼 远 超 吉 林 省 和 辽宁 省 。 由 此 推断 ,黑龙 江 省 水 稳 种 植 面 积 的 剧 增 
可 能 导致 了 该 省 农作物 受 旱灾 面积 与 SP 忆 指数 相关 性 的 降低 。 

(4) 由 于 土壤 湿度 数据 缺失 严重 采样 频率 过 小 且 易 受 耕 作 方 式 影响 等 原因 ,各 个 站 点 土壤 湿度 与 SPEI 
的 相关 系数 相对 偏 小 :在 w = 0.0 乞 的 显 彰 水平 上 (相关 系数 的 临界 值 为 0.17)，SPEIO1、SPEI03、SPEI06 和 
SPEI12 与 土壤 湿度 的 平均 相关 系数 分 别 为 0.43 .0.47 .0.48 和 0.44, 均 不 超过 0.5。 该 结果 可 以 从 赵 静 5) 的 研 
究 中 得 到 部 分 佐证 ; 闭 州 站 SPFEIol1 与 土壤 湿度 的 相关 系数 为 0.482, 稍 高 于 本 研究 的 0.468。 土 壤 湿度 和 
SPFI 指数 均 是 反映 干 湿 状 况 的 重要 指标 ,本 文 使 用 的 土壤 湿度 是 旬 值 数据 ,每 月 采样 3 次 ,而 SPEI 指数 则 是 
依据 逐日 的 气象 数据 转化 为 月 尺度 计算 而 来 ,其 精细 程度 远 超 土壤 湿度 数据 。 因 此 ,提高 二 者 相关 系数 的 关 
键 在 于 增强 十 二 湿度 数据 的 时 间 连 续 性 和 完整 性 ,如 提高 采样 频率 、 用 插 补 法 补 全 缺失 数据 等 。 


S 结论 


王 旱 指数 的 区 域 适用 性 是 地 区 干旱 精确 定量 化 的 前 提 和 基础 。 目 前 有 关 SPEI 指数 的 区 域 适用 性 研究 相 
对 缺乏 ,但 是 通过 我 们 统计 检验 以 及 相关 性 分 析 , 验 证 了 SPEI 指数 在 东北 地 区 干旱 研究 中 的 适用 性 。 本 文 主 
要 结论 如 下 : 

(1) 东 北 地 区 1、3、6 和 12 个 月 尺度 的 累积 水 分 亏 缺 量 序列 均 能 被 Log-logistic 分 布 函数 较 好 地 拟 合 , 基 
于 该 分 布 的 SPEI 指数 在 东北 地 区 的 应 用 具备 数学 统计 理论 基础 。 

(2)SPEI 指数 与 典型 干旱 事件 的 吻合 度 较 高 ,能 较为 准确 地 反映 干旱 事件 的 时 间 `. 地 点 和 强度 ,是 干旱 监 
测 和 预警 的 重要 手段 。 

(3)SPEI 指数 与 干旱 灾情 数据 呈 显 著 负 相关 。 干 旱 越 严 重 ,SPEI 指数 越 小 ,农作物 受 旱 灾 面 积 越 大 。 
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(4)SPEI 指数 与 土壤 水 分 资料 呈 显 著 正 相关 。SPEI 指数 能 指示 土壤 水 分 状况 ,SPEI 指数 越 大 ,土壤 湿 


度 越 大 。 


综 上 所 述 ,SPEI 指数 的 理论 假设 符合 东北 实际 情况 ,并且 能 被 历史 资料 和 事实 所 验证 ,说明 SPEI 指数 在 


东北 干旱 研究 中 的 适用 性 较 好 。 
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